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ABSTRACT 

For the first time, it has been shown experimentally that the Marangoni 
phenomenon is a primary mechanism for the movement of a gas bubble in a non- 
isothermal liquid in a low-gravity environment. In such a two-phase flow sys- 
tem, local variations in surface tension at the bubble surface are caused by a 
temperature gradient in the liquid. Shearing stresses thus generated at the 
bubble surface lead to convection in both media, as a result of which the bubble 
begins to move. A mathematical model consisting of the Navier-Stokes and 
thermal energy equations, together with the appropriate boundary conditions 
for both media, is presented. Parameter perturbation theory is used to solve 
this boundary value problem; the expansion parameter is the Marangoni num- 
ber. The zeroth, first- and second-order approximations for the velocity, 
temperature and pressure distributions in the liquid and in the bubble, and the 
deformation and terminal velocity of the bubble are determined. Experimental 
zero-gravity data for a nitrogen bubble in ethylene glycol, ethanol, and silicone 
oil subjected to a linear temperature gradient were obtained using the NASA 
Lewis zero-gravity drop tower. Comparison of the zeroth order analytical re- 
sults for the bubble terminal velocity showed good agreement with the experi- 
mental measurements. The first- and second-order solutions for the bubble 
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deformation and bubble terminal velocity are valid for liquids having Prandtl 
numbers on the order of one, but there is a lack of appropriate data to test the 
theory fully. 

INTRODUCTION 

The Marangoni ilow phenomenon, the flow generated by a gradient in sur- 
face tension, could be the driving flow mechanism in a low-gravity environ- 
ment. Because of the role of this flow phenomenon in materials processing in 
space, fluid storage and transfer operations in space, and boiling heat trans- 
fer as applied to propellant management in space, it is receiving increasing 
study. 

The phenomenon of surface tension induced flow requires a liquid-vapor 
interface and local variations in interfacial tension. In this two-phase flow 
system, local variations in the surface tension at the bubble Interface are 
caused by temperature variations in the liquid and gaseous media. According 
to the Marangoni flow phenomenon, liquid should flow around the bubble and the 
bubble should migrate toward the hotter region. However, it is difficult to 
study this flow phenomenon in normal gravity because of buoyancy effects. In 
normal gravity, buoyancy driven flows generally mask surface tension induced 
flows, that is, the magnitude of the buoyancy force is much greater than the 
thermophoretic force and buoyancy driven convection may further complicate 
the study. It is, therefore, virtually Impossible to carry out a quantitative 
experimental investigation for "pure" surface tension Induced flows in normal 
gravity. As a result, the experimental investigation of this flow phenomenon 
was conducted in the NASA Lewis Research Center Zero-Gravity Facility. 

The primary objective of this work was to observe the behavior of a single 
noncondensible bubble in a liquid under the Influence of a linear temperature 
gradient in a weightless environment, and then to compare these experimental 
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data with existing theory. A secondary objective was to extend the analytical 
model of this flow phenomenon beyond that which is presently available. 

Although the Marangoni flow phenomenon is a classical fluid physics prob- 
lem, only a limited amount of experimental and analytical work has been done 
in this area for the reasons mentioned above. Hershey in 1939 (ref. 1) conduc- 
ted an experiment in which it was shown that if a temperature gradient exists 
along a free liquid-vapor interface, a dynamic steady state is achieved. Block 
in 1956 (ref. 2) carried out experiments to prove that Benard (convection) cell 
motions are produced by variations in the surface tension due to thermal gra- 
dients. The first experimental work concerned with the study of surface ten- 
sion induced flow of a bubble in a heated liquid was done by Young, Goldstein, 
and Block in 1959 (ref. 3). Experimentally, they obtained a temperature gra- 
dient in a liquid sufficient to hold a bubble stationary in normal gravity by bal- 
ancing the gravity force with the thermophoretic force. Their experimental 
data was in fair agreement with their theoretical predictions. Hardy (ref. 4) 
repeated Young, et al.'s experiment in a more precise manner, and concluded 
that his data were in good agreement with the previous study. McGrew in 1968 
(ref. 5) conducted experiments to verify the theory he developed for computing 
the static thermophoretic force on a bubble. Papazian and Wilcox (ref. 6) were 
the first to conduct a Marangoni experiment in a low-gravity environment. 

Using a Space Processing Applications Research (SPAR) sounding rocket to 
achieve this type of environment, they were not successful in observing this 
flow phenomenon, primarily due to experimental difficulties. Theoretical analy- 
ses similar to that of Young, et al., were performed by Kuznetsov, et al. 

(ref. 7) and by Lyubln and Povitskly (ref. 8). Bratukhin's analysis (ref. 9) dif- 
fered from the previous investigators in that perturbation theory was used, 
where the dependent variables were expressed in a series expansion with the 
Marangoni number being the control parameter. 
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To date, relatively few quantitative experimental measurements of the 
Marangonl bubble motion flow phenomenon have been made to compare with the 
theory. The theory developed by the above investigators has been limited to 
the solution of the linearized Navier-Stokes equations (creeping flow approxi- 
mation) and the assumption of pure conduction which permitted the use of 
Laplace's equation, together with the appropriate boundary conditions on 
velocity, stress, and temperature. For a spherical bubble, the corrected ex- 
pression for the bubble terminal velocity as determined by Young, et al. is 

:kija 
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(See the symbol list for definitions of the quantities.) By setting the velocity 
equal to zero and neglecting the ratios of viscosities, densities, and thermal 
conductivities of the inner to outer media, the temperature gradient required 
to hold the bubble motionless can be found and is given by 


dz ;{ da/dT 

EXPERIMENTAL PROGRAM 


( 2 ) 


Initially, there were three key feasibility issues identified. The first had 
to do with the development of a bubble injection system to Introduce a single 
bubble into the liquid with minimal disturbance to the liquid and such that bub- 
ble size could be controlled. The second had to do with the development of a 
heating system to establish and maintain a linear temperature gradient in the 
liquid without causing any buoyancy driven convection. The third issue was an 
operational one, that is, the test run had to be completed within an S-hour time 
shift. Other concerns were the selection of test liquids, the determination of 
tank size, safety considerations, effect of bubble oscillations, selection of tem- 
perature instrumentation, visual obseiwation of bubble behavior, effect of sur- 
face contamination, drop-tower limitations, etc. 
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Based on the above mentioned key icusibility issues and other eonrerns, 
the experimental apparatus used to observe the Marangoni Bow phenomenon 
was designed around the Lewis Research Center's 155 m zero-gravity facility, 
shown in Figure 1. Data were obtained by allowing a capsule housing the ex- 
periment to free-fall for a time of 5.2 seconds. 

Experimental Apparatus 

A schematic of the experimental apparatus Is shown in Figure 2. A cy- 
lindrical test container, 12 cm In diameter and 12 cm high was fabricated from 
a transparent plastic material to permit visual observation of the bubble rising 
in the liquid. The container side walls were Insulated using a vacuum jacket 
fabricated from a transparent plastic material. The bubble injection system 
consisted of a supply of pressurized gas, a valving system consisting of a 
micrometer valve and a solenoid valve mounted in series, and a nozzle. A 
noncondensible bubble (nitrogen) was formed and introduced Into the liquid from 
a submerged orifice at the bottom of the test compartment. To establish and 
maintain the desired temperature gradient, a heater was employed at the top 
of the test container and a cooling bath at the bottom. Thermocouples attached 
to a rake, 2 cm apart and 3 cm from the tank centerline, were employed to 
monitor the liquid temperature before and during the test drop. Another ther- 
mocouple embedded In the heater plate was employed to monitor the plate tem- 
perature. The bubble movement was recorded by two high-speed cameras 
which were mounted 90° apart to observe the motion if the bubble path deviated 
from rectilinear. The elapsed free-fall time was obtained from two digital 
clocks In view of the cameras. A centimeter scale was mounted to the lower 
plate for scaling purposes. The bubble distortion resulting from the refractive 
index of the cylindrical tanks was accounted for when the data were analyzed. 

A photograph of the test apparatus Is shown In Figure 3. 
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Initial Tests - Feasibility Issues 

Preliminary zero-Rravit>’ tests were conducted to verify the bubble injec- 
tion system. The bubble was introduced into the liquid a fraction of a second 
after the drop was initiated. As a result of these tests, several experimental 
problems of Marangoni flow became apparent. First, a finite time for bubble 
oscillations to damp out was required, leaving fewer seconds of free-fall time 
to observe the Marangoni flow phenomenon. Second, the bubble had an appar- 
ent residual velocity even after 5.2 jcconds of free-fall time. To illustrate 
these problems, Figure 4 shows a representative plot of bubble displacement 
as a function of free-fall time for distilled water. The bubble residual velocitj' 
still persists after 3.2 seconds of free-fall time, but it appears to be decaying 
roughly exponentially. Also, note on the figure the free-fall time at which the 
bubble oscillations damped out as determined photographically. In other tests, 
methylene blue dye was introduced into the liquid to determine liquid disturb- 
ance as a result of bubble injection. These tests showed that the initial dis- 
turbance to the liquid decayed rapidly, and no gross motion of the liquid oc- 
curred. It was noted, however, that the dye affected the bubble formation and 
separation proc'esses, and it was impossible to form a single bubble. 

Besides the bubble injection system, another key feasibility issue identi- 
fied at the onset of this experin ent was whether or not the desired temperature 
gradient could be established and maintained in the liquid. The apparatus just 
described was designed with the intent of having the greatest flexibility for 
realizing this objective, without having any buoyancy driven convection present. 
Preliminary ground tests were conducted with distilled water to see if a linear 
temperature gradient could be obtained, aiul to determine the required heating 
time. Figure 5 shows a representative temperature profile after 260 minutes 
of heating. As can be seen, a linear temperature gradient is approached in the 
top half of the tank. The liquid was heated from the top to avoid Raleigh Insta- 
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blUties. Radial heat loss was minimized by the sidewall vacuum Jacket insu- 
lation. As a further check to insure there was no buoyar,v*y driven convection 
present, a separate thermocouple probe was used to measure the liquid tem- 
perature along the centerline iixis of the container. As the thermocouple was 
lowered Into the liquid through the standpipe, It took only a few seconds for the 
measurements to stabilize at each position. The results, when compared with 
the temperatures measured cm from the centerline, were nearly identical, 
indicating there was no buoyancy driven convection occurring within an annulus 
of fluid approximately 6 cm in diameter. 

Another key feasibility issue identified early in the program was the time 
required for a test. Although it was an operations problem rather than a tech- 
nical one, it was imperative that a test be completed with an 8-hour working 
shift. Time estimates of the various phases of the test program included 
ground preparation, transport of the expcriiv.ent vehicle, heating of the liquid, 
evacuating of the drop chamber, and recovery of the experiment. Summing up 
the maximum time estimates for each of the test stages, it was found that a 
complete test could be completed within an 8-hour time frame. 

Test Fluids 

Several factors were taken into consideration in the selection of the test 
fluids, the most important being the thermophysical properties, clarity, pre- 
sence of surfactants, and potential safety and health problems. Hie liquids 
chosen were distilled water, ethylene glycol, IX 200 silicone oil, and ethanol. 
Although distilled water is known to contain surfactants, it w.as chosen because 
of its high surface tension gradient and safety considerations, and was expected 
to provide experimental data outside the creeping flow regions. Kthylene glycol 
anil silicone oil were selected because any surface active agents present tend 
to stay in the bulk phase rather than migrate to the bubble surface. Kthanol was 
chosen because it is available as a reagent grade pure li(|uid, and experimental 
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data can be obtained outside the creeping How regime. A table of properties of 
these four test liquids is shown in Table 1, 

KXPERIMENTAL UKSULTS 

Meaningful data were obtained from about HO percent of the 38 tests per- 
formed in the 155 m drop tower. Complete results are presented by the author 
in reference 10. Single bubbles were generated in about half of the tests, but 
the tests when multiple bubbtes were formed (generally two) were also of value 
and arguments are presented to validate their usefulness. The importance of 
the no/.zle configuration is also discussed, as are bubble displacements, bubble 
si/e, shape, and oscillation measurements. The nozzle characteristics which 
were varied in the experiment included the nozzle height, that is, the distance 
measured from the bottom of the test container to the top of the nozzle, the in- 
si( e diameter, ami the nozzle material. 

The basic raw data collected from the experiments consisted of (1) taking 
bubble displacement measurements as a function of free-fall time from the 
films; (2) bubble shape measurements as a function of time; and (3) bubble os- 
cillations as a function of time. P'or accuracy of measurement and convenience, 
bubble displacement was measured from the orifice to the bottom of the bubble, 
while t!ie tin.c was measured from the initiation of free-fall. Bubble shapes 
were corrected for distortion using appropriate scale factors, while the bubble 
oscillations, after the initial vertical rise, were in the second mode, that is, 
the bubble oscillates between an oblate and prolate shape. 

Figure G shows the comparison of a typical isothermal run with a typical 
nonisothcrmal run for distilled water. As can be seen, the isothermal bubble 
displacement was continually decreasing with time, having a residual velocity 
of 0. 3 cm sec in the last second of free-fall. The nonoslthermal run shown is 
the one having the greatest displacement as a function of time with a bubble 
terminal \elocitv of 0.6 cm sec in the last second of free-fall when the tern- 


pcrature gradient was about C/cm. Slncv the expected value of xeloclty is 
an order of magnitude greater than this, there is very little discernible differ - 
ence between the velocities obtained from these curves, and it appears that 
the Marangcni flow phenomenon was not observable using distilled water. In 
addition, it should be noU'd thnt a single bubble was not generated In any of 
these runs. Marangoni flow was not observed In distilled water. This may 
have been because of surfaetants concentrated at the bubble surface or because 
of some molecular phenomenon, perhaps caused by the polar nature of the 
water molecule, that occurred at the bubble interface. Either of these factors 
could eliminate or change the surface tension gradient. 

For the other three test liquids, ethylene glycol, DC 200 silicone oil, and 
ethanol, Marangoni flow was observed. Figure 7 is representative of the ex- 
perimental data obtained for the three test liquids. In this figure, data are pre- 
sented for ethylene glycol. Bubble displacement is plotted as a function of 
free-fall time for four temperature gradients ranging from 1.0° to 4.0° C/cm. 
For comparative purposes, results for an isothermal run are also shown. The 
results show that the bubble displacement for the nunisothermal cases becomes 
nearly linear with time, indicating that the bubble reached a terminal velocity. 
The bubble displacement for the isothermal case becomes relatively flat, hav- 
ing a residual velocity of less than 0.05 cm /sec over the lust second of free- 
fall time. Bubble diameter for each of these test runs was approximately 0.8 
cm. The bubble shape, after the oscillations damped out In I;bs than 0.2 sec- 
onds, was spherical. Some of the test data when multiple bubbles were formed 
was used. Although the bubbles are follou'ing mutually perpendicular paths, 
arguments supported by experimental data indicate the usefulness of the data 
dep«‘nding on the relative sizes of the bubbles .and their separation distance, 
as shown by the author (ref. 10). Some of the more important multiple-bubble 
data results are now breifly discussed. 
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In this Htu(ly, the prim'iplc Interest was in observinji the classical-fluid 

I 

physics Mnranftoni flow phenomenon for u single bubble in a r.ero-gravity envi- 
ronment. The experimental bubble terminal velocity and bubble shape data 
were to be compared with existing theory, and an extension to the theory devel- 
oped as part of this dissertation. However, in the actual rase of materials 
processing in space, multiple bubbles will likely be present. Consequently, 
future investigations of the Marangoni bubble phenomenon should include 
studies of multiple bubble interaction effects as well as a continuation of single 
bubble experimentation. As determined from the experimental data, important 
multiple bubble interaction factors include: the relative bubble sizes, whether 
a bubble is leading or trailing another bubble, bubble separation distance, and 
relative bubble paths, that is, whether the paths arc perpendicular, parallel, 
or at some other angle. Some of the conclusions reached from the limited 
amount of experimental data obtained Indicate that bubble interaction effects is 
a very worthwhile area to pursue. For example, for two bubbles, it was de- 
termined that the leading bubble wake does affect the temperature distribution 
field around the trailing bubble, especially when the paths are perpendicular. 

In this case, if the leading bubble is the larger one, its terminal velocity does 
not appear to be affected, whereas the trailing bubble velocity is greatly affec- 
ted and, in some cases, the trailing bubble comes to rest. On the other hand, 
when the trailing bubble is the larger one, it rises faster, indicating the wake 
effect of the leading bubble is minimal on the trailing bubble, and the bubbles 
should coalesce. In the cases where the bubble paths were not perpendicular, 
there appeared to be little bubble interaction effects even when the leading 
bubble was larger than the trailing bubble. Of course, the bubble sizes as well 
as bubble separation distance also affect the bubble terminal velocities. 

In addition to the bubbic-bubble interaction phenomenon, another interes- 
ting phenomenon which was obsen’ed in practically all the isothermal zero- 
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gravitv rcfltfi was that, regardless of the test liquid used, the bubbles after 5.0 
seconds of zero-gravity time still had some residual velocity. Longer zero- 
gravit>' times are required to determine how iong it tvould take for the bubble 
to approach zero velocity since the only force acting on the buU'le interface 
is the drag force. The conclusion reached from the experimental results is 
that the residual velocity, after 5.0 seconds of free-fall time, is inversely 
dependent on the viscosity of the liquid. For all test liquids used, the bubble 
velocity decay rate appears to be exponential. 

Also, it was observed that the bubble size did not change during the 5.0 
seconds of free-fall time, which indicates that there was no observable bubble 
expansion during this time due to the heating of the bubble. To study the heat- 
ing effect on bubble expansion would require much longer zero-gravity times 
than can be obtained in a ground-based research facility and, therefore, this 
would have to he conducted in a space environment. 

Some representative data for thi' three test liquids are shown in Fig- 
ure H. The linear temperature gradient for the silicone oil was on the order 
of 2° C/cm, whereas for ethylene glycol and ethanol, it was on the order of 
4° C/cm. Ihc circles represent the relative bubble sizes, which for sili- 
cone oil and ethylene glycol were on the order of 0.8 cm, while for •ithanol 
the size was approximately 0.6 cm. A teflon tip nozzle was used for the 
ethanol as opposed to a stainless steel tip nozzle for the other two fluids. 

Note from the figure the nearly linear displacement as a funi tion of time for 
all three lii|uids, indicating in each case that the bubble reached terminal 
velociiy . 

As an illustration of tin* role of the nozzle configuration in the formation 
of u single bubble, thi' following discussion is presented for ethylene glycol 
as the test fluid. Initially, a nozzle height of 9 cm was chosen, but multiple 
bubbles were always formed by a stainless steel nozzle. In addition, it was 
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observed that the proximity of the top plate affected the bubble diapiacement. 

As a eonaequenc'c of having these undesirable conditions, a nozzle height of 
6 cm was chosen, which again resulted in multiple bubbles. A teflon tip was 
placed over the stainless steel nozzle with the same diameter in an effort to 
form a single bubble. This solution was partially successful in that a single 
bubble was formed, but not until late into the drop. The late bubble forma- 
tion was attributed to the relative difference in the contact angle betw’cen the 
teflon and t^e stainless steel. A final attempt to produce a single bubble was 
made using a smaller diameter orifice. Decreasing the orifice inside diameter 
from 0.066 to 0.043 cm and a height of 6 cm enabled three single bubble exper- 
iments to be conducted. 

The experimentni data obtained in th>s study are correlated in terms of 
two dimensionless parameters, the Keynolds number and the Marangoni num- 
ber. The Keynolds number, a ratio of inertial to viaeous forces, is given by 

P,ar 

Kc = -i— (3) 

and the Marangoni number is given by 
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whi're I' > bubble terminal velocity, a bubble radius, liquid viscosity, 

- liquid density, Ao AT = surface tension gradient, and AT/Az tempera- 
ture gradient. 

ANALYSIS 


In developing the matiH'mutieuI analysis for the problem, the following 
assumptions are m:ide. The problem is formulated as being sU'ady state with 
uniform streaming flow ami a linear temperature gradient at infinity, .\xial 
symmetry exists in spherical coordinates and the bubble is allowed to lieviate 
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slightly from spherical. The two media, gas and liquid, arc assumed to be 
Immiscible, and the liquid is Newtonian. The force of gravity is neglected 
and the thermophysical properties of the two media are constant, that is, den- 
sity, specific heat, thermal conductivity, and viscosity, even though the temper- 
ature is -irying. The one quantity allowed to vary is the surface tension a, 
which, for most liquids, depends linearly on temperature. Mathematically, the 
surface tension variation can be expressed as 

cr-o +^AT (5) 

° dT 

where (t^ is a reference surface tension at temperature T^ and the surface 
tension gradient, do/dT, is a constant value. 

The governing equations for this problem are the Navier-Stokes equations, 
the energy equations and the boundary conditions.. Boundary conditions are 
applied at the bubble interface to express (a) the imjjenetrability (no flow) con- 
dition of the normal velocity; (b) the continuity of the tangential velocities; 

(c) the continuity of normal stress; (d) the continuity of tangential stress, in- 
cluding that due to the therma’. variation in surface tension; (e) the continuity of 
temperature and heat flux; and (f) the finiteness of the variables at the limits of 
zero and infinity for the rad'al variable. 

The next step in the analysis was to noiidimensionalize these expressions, 
using a minimum parametric representation. The three basic scale factors de- 
termined from this method were the stream function, pressure, and velocity 
given by 
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Substituting the dimensionless ratios with the scale factors P^, 
and into the governing equations and boundary conditions, the following 
dimensionless equations result: 
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Normal stress condition at the interface 
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Tangential stress rondition at the inUTlace 
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Temperature continuity at the interface 
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Condition ut I'initcness as radius appraachcs zero 

'l<„ -► finite 

Condition at uniform How as radius approaches infinity 
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Boundedness conditions on temperatures 
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The dimensionless groupings appearing in the above equations are: 
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Thev are defined as: 
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whieh is ealled the Marangoni number, Pr is the Prandtl 
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number, — , — , — are the dimensionless ratios of density, viscosity, 
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thermal conductivity and specific heat, and, finally. 
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All of the above dimensionless groupings are calculable, based on known fluid 
thermophvsical properties, a temperature gradient, and a characteristic length. 
It Is noteworthy to mention that once the lluids are chosen, there arc only t\\’o 
control parameters remaining, namely the temperature gradient, dT^/dz, and 
the characteristic length, usually the bubble radius, both of which appear as part 
of the Marangoni number. The general ecjuations governing the motion of a 
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deformable bubble, including the boundary conditions in a nonisothermal liquid, 
have now been derived in dimensionless form. The next step in the analysts is 
to obtain analytical solutions . 

Parameter Perturbation Theory 

Parameter perturbation theory was used to solve the governing equations 
together with the appropriate boundary conditions. In this theory we seek a so- 
lution to the boundary value problem in terms of a control parameter, that is, 
the dependent variables are expanded in terms of this control parameter. Basic- 
ically, what is done to determine approximate solutions to differential equations, 
whether they be governing equations or boundary conditions, is to assume ex- 
pansions, substitute them into the equations, and then solve the equations. An 
important class of such asymptotic sequences is represented as a collection of 
integer powers of the perturbation parameter, in this case, the Marangoni num- 
ber. The general expression for a variable f(r, Ma) is given by 

00 

f(r, Ma)=^ ajj(r)Ma” (20) 

n=0 

as the Ma approaches zero, where a^(r) is an expansion variable. This is a 
straightforward expansion of the Poincare* type and is used in this analysis. 

Care must be taken to insure that the expansion is uniformly valid. As a re- 
sult of making these substitutions, sets of equations for the zeroth-, first-, and 
second-order equations for both the governing equations and the boundary condi- 
tions arc obtained. Solutions for these sets of equations are obtained. Young, 
et al. (ref. :i) obtained the zeroth-ordcr solution, and Bratukhin (ref. 9) pre- 
sented the first-order approximation. Second-order solutions were obtained in 
reference 10 for the stream functions, temperature, and pressure distributions 
in the inner and outer media, and expressions for the bubble terminal velocity 
and the bubble deformation. In dimensionless form, the bubble deformation and 
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the terminal velocity (in terms of Reynolds number) are 

R= 1 _ _2_ Pr\ I Ma^ 

Uq ,128 144 \ 2 j 
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- :t cos 0 
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( 21 ) 


and 


Re 
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Vl28 144 1 


9 
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+ _L pr - i£L Pr2 


320 


8,640 


Ma' 


( 22 ) 


In this case the bubble radius is a function of the Prandtl number and the dimen- 
sionless parameter for both the first- and second-order approximations. 

On the other hand, the terminal velocity is affected by the zeroth- and the 
second-order approximations, that is, there is no contribution to the terminal 
velocity resulting from the first-order approximation. Here too, the second- 
order effect is a function of the Prandtl number and the parameter 

Surface tension variation must be included in this analysis. It is important 
to note that the Navier-Stokes equations and the energy equations are thereby 
coupled through the normal and shear stress boundary conditions which contain 
surface tension terms, and the surface tension is a function of the liquid tem- 
perature. 

It was shown from the analysis that no nonuniformities result from the sub- 
stitution of lowcr-ordcr solutions into the governing equations and boundary con- 
ditions to obtain higher-order solutions. Thus, the dependent variables can be 
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expanded In a straightforward manner In terms of Integer powers of the 
Marangoni number. Generally, the nonuniformity condition which arises In 
isothermal bubble motion in normal gravity results because the uniform veloc- 
ity condition as | r| approaches infinity cannot be satisfied in the first ap- 
proximation. This is commonly referred to as Whitehead's paradox. In the 
present study, this did not occur in the first- and second-order approxima- 
tions. However, based on the results of the second-order approximation, it 
is suggested by the author (ref. 10) that higher-order approximations, diird 
and fourth order, may result in the development of a nonuniformity condition, 
that is, the boundary condition of uniform velocity at Infinity will not be satis- 
fied. This conclusion Is based on the presence of the Skokeslet term, which 
is the r term, in the second-order approximation of the stream function. 

This condition on the straightforward expansion may not appear until the 
fourth-order approximation, which should contain the next approximation to 
the bubble terminal velocity. At some higher order, the boundary condition 
at infinity cannot be satisfied. At that point another solution technique would 
have to be used, such as a matched asymptotic expansion, integral, or num- 
erical method. 

COMPARISON OF EXPERIMENTAL AND THEORETICAL RESULTS 
Experimental and theoretical results for ethylene glycol and silicone oil 
are plotted in Figure 9. For the experimental data, the Marangoni number is 
plotted against the Reynolds number for single bubbles. Theoretical results 
are plotted for both the zeroth-order and the higher-order solutions . The 
zeroth-order solution (creeping flow), given by the first term in equation (22), 
is shown as the solid line. As was determined earlier, the higher-order so- 
lution, which includes the second-order approximation given by the second 
term in equation (22), has a very limited range of applicability, even for the 
limiting case where the Prnndtl number approaches zero and the dimensionless 
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paruinetcr approaches infinity. ' 11)0 dashed line (theoretical result) shown 
in Fii^rc 9 illustrates this case Kraphicallyt it deviates from the experimental 
results very sharply for Reynolds numbers greater than 1.5. What is, how- 
ever, most noteworthy from this figure In comparing experimental and theo- 
retical results, is that the experimental data generally match the zeroth-order 
solution, even for the relatively high Marangoni numbers. In other words, for 
Marangonl numbers greater than 0.1, the theory is not expected to be appli- 
cable, but it apparently is. From the first-order approximation, it is predic- 
ted that, Independent of the Prandtl number, a standing eddy should begin to 
form behind a bubble when the Marangoni number is approximately equal to 
16/3. As the standing eddy begins to develop, it should result in changing the 
surface tension gradient over the rearward portion of the bubble interface. If 
this occurs, the experimental results should begin to deviate from the theo- 
retical results at Marangoni numbers greater than 16/3. Comparing this ana- 
lytical result with the experimental data in Figure 9, the data show that there 
is \ery little deviation from the zeroth-order solution up to Marangoni numbers 
of about eight. Substitution of typical values for Prandtl number and for 
these test liquids into the expression for first- and second-order bubble de- 
formation predicts that the bubble should change shape for intermediate values 
of Marangoni numbers. However, in every test, the bubble shape, after the 
initial oscillations damped out, remained spherical. This suggests that the 
equation for bubble deformation .s valid for liquids having Prandtl numbers on 
the order of one. The number of data points is rather limited and before any 
definitive conclusions can be drawn from the comparison of the experimental 
and theoretical results, additional data points are needed. In particular, if 
comparisons are to be made at very low Marangoni numbers using relatively 
high viscosity test liquids, much longer zero-gravity times than the 5.2 sec- 
onds obtained in the drop tower are required. In faet, it would also be 
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desirable to have longer zero-gravity times over which to take data for the 
higher Marangonl numbers. 

A compilation of all of the experimental data, Including the multiple- 
bubble as well as the single-bubble data, is presented in Figure 10 for the 
three test liquids in which Marangonl flow was observed. Jata points for mul- 
tiple bubbles are included on this figure since arguments were made earlier to 
validate their usefulness. For the multiple-bubble data points, the subscripts 
L and T arc used, where the subscript L indicates a leading bubble and the 
subscript T indicates a trailing bubble. For comparison, the zeroth-order 
approximation is shown as the solid line in the figure. The dashed line through 
the data points indicates that at the higher Marangonl numbers, the experimental 
data does deviate from the theory. One possible explanation for the deviation is 
that perhaps a standing eddy is being formed behind the bubble, as predicted 
from the first-order solution mentioned previously. As the standing eddy 
grows, the surface tension gradient is affected more and more, thus causing 
the bubble terminal velocltj' to decrease, which correspondingly reduces the 
Reynolds number. As was determined earlier, the bubble shape does not affect 
the bubble terminal velocity. The bubble remains spherical, in contradiction 
to the theory. The data shown in this figure prove the existence of the 
Marangonl flow phenomenon, and shottld serve as a basis for a continuing re- 
search effort in this area, both experimentally and analytically. 

CONCLUSIONS AND IlKCOMMENDATIONS 

The principle conclusion of this research is that the Marangonl bubble mo- 
tion flow phenomenon does occur and can be a driving flow mechanism in a 
zero-gravity environment. For the first time, it has been shown experimentally 
that local surface tension variations at the bubble Interface caused by a tempera- 
ture gradient in the liquid cause the bubble to migrate from the cold to the hot 
region and thus the flow phenomenon is real. This flow phenomenon was ob- 




served in three out of the four test liquids used. The range of Rejnolds nuin* 
bers was from about 0.25 to 00, with corresponding Marangoni numbers 
ranging from about 0.50 to 600. This flow phenomenon was not observed in 
distilled water. 

Comparison of the experimental with the theoretical results were made on 
the basis of bubble terminal velocity. Experimental data could be correlated 
in terms of two parameters, namely the Reynolds number and the Marangoni 
number. In the zeroth-order approxi nately solution, the Reynolds number 
was proportional to the Marangoni number, and It was proved analyticolly 
that the bubble must remain spherical. Although the theory In this case would 
be expected to be valid only for Marangoni numbers much less than one, the 
theory was in good agreement with experimental data for Marangoni numbers 
up to about eight, at which point the data begins to deviate from theory. Ex- 
perimental data also corroborated the analytical result that the bubble remains 
spherical. Not only was this the case for small Marangoni numbers, but for 
Marangoni numbers approaching 600 as well, which contradicts higher-order 
theory. 

A dimensional analysis was performed on the governing equations which 
include the Navier-Stokes equations and the energy equations together with the 
appropriate bounoary conditions for the inner and outer media. The following 
dimensionless ratios and groupings evolved: viscosity, density, thermal con- 
ductivity, and specific heat ratios of the inner and outer media, aiiu the dimen- 
sionless parameters, the Marangoni number and the Prandtl number, expressed 
in terms of the outer medium fluid properties. Two other dimensionless group- 
ings are the uniform velocity of the liquid at infinity and a grouping expressed as 


. Once the Inner and outer media are choeeu, Jiere are only two 



control parameters remaining, namely, the temperature gradient and the 
characteristic length, which in this case is the bubble radius. Both of these 
parameters appear in the Marangonl number. For this reason, it was conclu- 
ded that the control parameter that should be uned in parameter perturbation 
theory was the Marangoni number. In other words, all the dependent variables 
are expanded in terms of the Marangoni number. 

In the first-order approximation, the bubble deformation depended not only 
on the Marangoni number, but also on the Prandtl number and the dimeni>ton- 
less parameter o^. Bubble deformation, however, was not experimentally ob- 
serA'ed, and this disagreement is due to the theory only being valid for the lim- 
iting case of Prandtl numbers less than one and values of the dimensionless 
parameter approaching infinity. The bubble terminal velocity was the same 
as in the zeroth-order approxination. 

From the first-order approximation, independent of the Prsndtl number, 
it was also found the', a standing eddy should begin to form behind a bubble 
when the Marangoni number is approximately equal to 16/3. r'vs the eddy de- 
velops, it should result in changing the surface tension gradient over the rear- 
ward portion of the bubble. If the eddy does oceur, the expvvimenial data 
should begin to deviate from the theoretical results at a Marangoni number 
somewhat greater than 16/3. The data indeed begins to deviatu from theory 
in the Marangoni range of 6 to 10, which is consistent with thftor>'. 

From the second-order approximate solution, it was dcterinincii that there 
was a second-order effect on both the bubble shape and its terminal velocity. 

The solutions for predicting second-order effects on bubble shape and bubble 
terminal velocity are valid if the Prandtl number is on the order of one and 


the (liinenHlonles parameter cr^ is Urge. Unfortunately the data obtained 
were for liquids having Prundtl nunibern much greater than one, and, therefore, 
die theory cannot be (11ri*ctly eomparci) with the data. It la, however, worth 
mentioning that one area where the areoml-order approxlmatlona for bubble 
terndr.al velocity, or Jteynoldp number, and bubble deformation may have 
aomt i'.pplication is for liquid metrla. In the caae of liquid metala, the Prandtl 
numbera arc mucii Icita than one. Marangoni bubble motion experimenta would 
have to be performed In a apace environment, auch aa in Spacelab, where long 
zero-gravity times arc available. 

The experimental data obtained in the preaent study prove the existence of 
the Marangoni flow phenomenon for bubble motion, and should serve ap a basis 
for a continuing research. The following arc areas where continuing research 
efforts could profitably be directed: 

1. It would be desirable to obtain data points for several test liquids in or- 
der to evaluate the effect of Prandtl number and the dimensionless parameter, 
0 ^, on the bubble terminal velocity and shape, as suggested from higher-order 
approximate solutions. Although the experiment does not agree with the higher- 
order approximations, qualitatively, the bubble terminal velocity and bubble 
shape may depend on a functional relationship of these two dimensionless pa- 
rameters. In other words, different curves of Reynolds .'-umber versus 
Marangoni -lumber may be obtained us and Prandtl number are varied. 

2. It would also be desirable to obtain data for u given Marangoni number 
by varying the two control variables, namely, the temperature gradient and the 
bubble radius, to sec if the bubble terminal velocity, or correspondingly, the 
Reynolds number is affected. For example, us the bubble radius is decreased 
below 1 centimeter, the temperature gradient must be increased in proportion 
to the radius squared in order to keep the Marangoni number a constant. Of 
course, this means that much higher temperature gradients would be required 
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than thoae obtained as part of this experiment, and this would require much 
longer heating times. 

3. Additional data should be obtsined at very low Marangonl numbers (the 
creeping flow regime), so that theory can be comp^arc^l with experimental data. 
This could be accomplished by using relstlvely low vlscoiUty fluids, such as 
silicone oil, and would require much longer zero-gravity times than the 5.2 
seconds obtained in the drop tower. It would also be desirable to have longer 
zero-gravity times at the higher Marangonl numbers because of uncertainties 
about whether the bubble formation and separation processes disturb the tem- 
perature field. 

4. Data should be obtained for bubble expansion resulting from heating by 
the llqulo. Much. longer zero-gravity times than those obtained In the drop 
tower are required. In all the nonisothcrmal experiments conducted in the 
drop tou'cr, there was no discernible change in bubble size. 

o. To extend the perturbation analysis to higher-order approximations 
does not appesr to be profitable at this time until low Prandtl number dsta are 
obtained. Perhaps better higher-order approximations could be achieved 
using parameter perturbation theory to match the Intermediate Marangonl 
number experimental data by expanding the dependent variables about a con- 
trol or expansion parameter other than the Marangonl number, for example, 
the Prandtl number times the Msrangonl number. 

6. Aiwther loglcsl step in extending the analysis Is to solve the governing 
equations, together with the boundsry conditions presented In this paper, nu- 
merically. If this Is achieved, the final step would be to allow the thermo- 
physical properties to be temperature dependent. 
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SYMBOLS 

a spherical bubble radius 

a^(r, 0) expansion variable 

C specific heat 

P 

D scalar differential operator 

f function 

g gravitational acceleration 

H mean radius of curvature 

h thermal conductivity 

K thermal conductivity 

Ma Marangoni number 

n integer 

Pr Prandtl number 

R deformed bubble radius 

Re Reynolds number 

r radius in spherical coordinates 

T temperature 

dT^/dz temperature gradient at infinity 
U bubble terminal velocity 

z |X)lar axis 

0 polar angle in spherical coordinates 

P 


viscosity 
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p density 

ff reference surface tension 

o 

dimensionless group 
da/dT surface tension gradient 

♦ stream function 

Subscripts: 

1, 2 outer and inner media 

o scale factor 

Superscript: 

( )' derivative 
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Figure 1. - NASA zero-gravity faciiity. 





Figure 3, - Zero-gravity experimental apparatus. 
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